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Double real radiation corrections to gluon scattering at NNLO 
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We use the antenna subtraction method to isolate the double real radiation infrared singularities present in 
the six-gluon tree-level process at next-to-next-to-leading order. We show numerically that the subtraction term 
correctly approximates the matrix elements in the various single and double unresolved configurations. 
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1. Introduction 

In proton-proton collisions, the factorised form 
of the inclusive cross section is given by, 
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(1) 



where d&ij is the parton-level scattering cross sec- 
tion for parton i to scatter off parton j normalised 
to the hadron-hadron flu^H and the sum runs over 
the possible parton types i and j. The probabil- 
ity of finding a parton of type i in the proton 
carrying a momentum fraction £ is described by 
the parton distribution function fi(£)d£. By ap- 
plying suitable cuts, one can study more exclu- 
sive observables such as the transverse momen- 
tum distribution or rapidity distributions of the 
hard objects (jets or vector bosons, higgs bosons 
or other new particles) produced in the hard scat- 
tering. The leading-order (LO) prediction is a 
useful guide to the rough size of the cross sec- 
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2 The partonic cross section normalised to the parton- 
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of £i and £2 into d&ij . 



tion, but is usually subject to large uncertainties 
from the dependence on the unphysical renormal- 
isation and factorisation scales, as well as possi- 
ble mismatches between the (theoretical) parton- 
level and the (experimental) hadron-level. 

The theoretical prediction may be improved by 
including higher order perturbative predictions 
which have the effect of (a) reducing the renor- 
malisation/factorisation scale dependence and 
(b) improving the matching of the parton level 
event topology with the experimentally observed 
hadronic final state pQ. The partonic cross sec- 
tion d&ij has the perturbative expansion 
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where the next-to-leading order (NLO) and next- 
to-next-to-leading order (NNLO) strong correc- 
tions are identified. For many processes, knowl- 
edge of the NLO correction is sufficient. However, 
for the main 2 — > 1 or 2 — > 2 scattering pro- 
cesses such as Drell-Yan, Higgs production, di-jet 
production, vector-boson plus jet, vector-boson 
pair production or heavy quark pair production, 
the NLO corrections still have a large theoret- 
ical uncertainty and it is necessary to include 
the NNLO perturbative corrections. In addition 
to reducing the renormalisation and factorisation 
scale dependence there is an improved matching 
of the parton level theoretical jet algorithm with 
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the hadron level experimental jet algorithm be- 
cause the jet structure can be modeled by the 
presence of a third parton. In this talk, we are 
mainly concerned with the NNLO corrections to 
di-jet production where the resulting theoretical 
uncertainty is estimated to be at the few per-cent 
level BQ. 



2. Antenna subtraction 

Any calculation of these higher-order correc- 
tions requires a systematic procedure for extract- 
ing the infrared singularities that arise when one 
or more final state particles become soft and/or 
collinear. These singularities are present in the 
real radiation contribution at next-to-leading or- 
der (NLO), and in double real radiation and 
mixed real-virtual contributions at next-to-next- 
to-leading order (NNLO). 

There have been several approaches to build a 
general subtraction scheme at NNLO [213141516171 
I8|9)10|ll)12j . We will follow the NNLO antenna 
subtraction method which was derived in [2] for 
processes involving only (massless) final state par- 
tons. This formalism has been applied in the com- 
putation of NNLO corrections to three-jet pro- 
duction in electron-positron annihilation [131141 
115116] and related event shapes [17118119120(21] . 
It has also been extended at NNLO to include one 
hadron in the initial state relevant for electron- 
proton scattering 22 23] while in refs. [24125] the 
extension of the antenna subtraction method to 
include two hadrons in the initial state is dis- 
cussed. 

In the antenna subtraction method the sub- 
traction term is constructed from products of 
antenna functions with reduced matrix elements 
(with fewer final state partons than the origi- 
nal matrix element), and integrated over a phase 
space which is factorised into an antenna phase 
space (involving all unresolved partons and the 
two radiators) multiplied with a reduced phase 
space (where the momenta of radiators and un- 
resolved radiation are replaced by two redefined 
momenta) . The full subtraction term is obtained 
by summing over all antennae required for the 
problem under consideration. In the most general 



case (two partons in the initial state, and two or 
more hard partons in the final state), this sum 
includes final-final, initial-final and initial-initial 
antennae. 

The relevant antenna is determined by both the 
external state and the pair of hard partons it col- 
lapses to. In general we denote the antenna func- 
tion as X. For antennae that collapse onto a hard 
quark-antiquark pair, X = A for qgq. Similarly, 
for quark-gluon antenna, we have X = D for qgg 
and X = E for qq'q' final states. Finally, we 
characterise the gluon-gluon antennae as X = F 
for ggg, X — G for gqq final states. At NNLO we 
will need four-particle antennae involving two un- 
resolved partons and one-loop three-particle an- 
tennae. 

In all cases the antenna functions are de- 
rived from physical matrix elements: the 
quark-antiquark antenna functions from 7* — > 
qq + (partons) [27J, the quark-gluon antenna 
functions from x ~ * 9 + (partons) [28) and the 
gluon-gluon antenna functions from H — > (par- 
tons) [29] . 

A key element in any subtraction method, is 
the ability to add back the subtraction term inte- 
grated over the unresolved phase space. In the an- 
tenna approach, this integration needs to be per- 
formed once for each antenna. Although no prob- 
lem in principle has been identified, at present, 
not all of the necessary integrals have been com- 
pleted. The current state of play is summarized 
in Table 1. 

3. Double real radiation gluonic contribu- 
tions 

We now consider the six-gluon contribution 
to the NNLO di-jet cross section. At leading 
colouiH, the double real radiation contribution is 
given by, 
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3 To simplify the discussion, we systematically drop the 
sub-leading colour terms. 
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Table 1 



Integrated antenna functions 





final-final 


final-initial 


initial-initial 
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sm 


(in progress) [H] 
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(in progress) 


(in progress) 
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(in progress) 



where the normalisation factor AT includes all 
QCD-independent factors as well as the LO de- 
pendence on the renormalised QCD coupling con- 
stant a s . S n /Z n is the group of non-cyclic per- 
mutations of n symbols and Ag(er(l), . . . , er(6)) 
is the square of a colour ordered amplitude, p\ 
and P2 are the momenta of the ingoing gluons, 
while ps , . . . , fig are the momenta of the four final 
state gluons and where the phase space integra- 
tion d$4 is over the 4-parton final state provided 
that precisely 2-jets are observed by the jet algo- 
rithm J2 . Here singularities occur in the phase 
space regions corresponding to two gluons becom- 
ing simultaneously soft and/or collinear and the 
subtraction term must successfully regularize all 
double unresolved and single unresolved singular- 
ities. 

The antenna subtraction term dafj NLO for the 
six-gluon process was derived in Ref. [2 5) . 

4. Numerical checks 

In this section we will test how well the sub- 
traction term d&pf NLO approaches the double real 
contribution d&^ NLO in the single and double 
unresolved regions of phase space so that their 
difference can be integrated numerically over the 
four-dimensional phase space. We will do this 
numerically by generating a series of phase space 
points that approach a given double or single un- 
resolved limit. For each generated point we com- 
pute the ratio, 



dF — 



(4) 



where da^ NLO is the matrix element squared 
given in equation ([3]) and dafj NLO is the sub- 
traction term of Ref. 125 . The ratio of the ma- 



trix element and the subtraction term should ap- 
proach unity as we get closer to any singularity 
and the quality of this convergence represents the 
first sanity check on the NNLO antenna subtrac- 
tion method for the six-gluon process. 

For each unresolved configuration, we will de- 
fine a variable that controls how we approach the 
singularity subject to the requirement that there 
are at least two jets in the final state with pt >50 
GeV. The centre-of-mass energy y/s is fixed to be 
1000 GeV. All of the single and double unresolved 
regions of phase space were analysed separately. 
In each case, we generated 10000 random phase 
space points and computed the average value of R 
together with the standard deviation. The results 
of the numerical study for all limits (apart from 
the single collinear limit we shall discuss below) 
are summarized in Table 2. 

The numerical results show that the combina- 
tion of the antennae in the subtraction term cor- 
rectly describes the infrared singularity structure 
of the matrix element. All of the double unre- 
solved and single soft singularities present in the 
matrix elements are cancelled in a point-by point 
manner. 



Finally we generate points corresponding to the 
final and initial state single collinear regions of 
the phase space. For the final-final collinear sin- 
gularity, fig. [Ha) shows the ^-distribution for a 
range of values of x = Sy/s. Similarly in the 
initial-final collinear limit, we show the distribu- 
tions of R for small values of x = su/s in fig.[2ja). 
We see that the subtraction term, which is based 
on azimuthally averaged antennae functions does 
not accurately describe the azimuthal correlations 
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Table 2 

Results for the average and standard deviation of R for 10000 phase space points close to singular limits. 





small parameter 


< R > 


(J 


double soft 


(s - 8ij)/a = 10" 6 


0.9999994 


4 x 10" 5 


triple collinear (final) 


Sijk/s = 10~ 9 


1.0000004 


4.2 x 10~ 5 


triple collinear (initial) 


sijk/s = -10- 10 


0.99954 


0.04 


soft/collinear (final) 


(s - s ijk )/s = 10- 6 , Sij /s= 10- 6 


0.99999993 


0.0001 


soft/collinear (initial) 


(s - s ijk )/s = 10- 7 , su/s = -10- 7 


0.99999998 


1.6 x 10~ 7 


double collinear (final-final) 


8ij /a= 10" 8 , s kl /s= 10- 8 


0.9999995 


0.00037 


double collinear (final-initial) 


s jfc / s = 10- 10 , = -10- 10 


1.00012 


0.018 


double collinear (initial-initial) 


s u /s =-10- 10 , s 2j /s = -10- 10 


1.00001 


0.004 


single soft 


- s tjk )/s = 10- 7 


0.999998 


1.9 x 10" 5 



present in the matrix elements and antenna func- 
tions when a gluon splits into two collinear glu- 
ons. This is the reason why the distributions in 
figs. [U(a) and [2ja) have such a broad shape. It 
is clear that as we approach the collinear limits 
x 0, the azimuthal terms are not suppressed 
and the subtraction term is not, point by point, 
an adequate representation of the matrix element. 

Nevertheless, the azimuthal terms coming from 
the single collinear limits do vanish after an az- 
imuthal integration over the unresolved phase 
space. Here we are performing a point-by-point 
analysis on the integrand defined by the ma- 
trix element squared and the subtraction term. 
One possible strategy is to introduce a correc- 
tion term Q F o(i,j, z, kx) to the F® four-gluon 
antenna functions which reconstructs the angu- 
lar terms [25]. Subtracting Q F o(i, j, z, k±) from 
the final-final and, initial-final and initial-initial 
configurations (by crossing momenta to the initial 
state) produces a subtraction term that is locally 
free of angular terms. 

With this azimuthally modified subtraction 
term, we recompute the distributions in figs.[ljb) 
and EJb). For 10000 final state single collinear 
phase space points and x — 10~ 12 we obtained 
an average R = 0.99994 and a standard deviation 
of (j — 0.015. We repeated the same analysis for 
the initial state collinear configuration for x = 
— 10~ 12 and obtained an average of R = 1.00007 
and a standard deviation of a = 0.012. For both 
cases, the distributions now peak around R = 1 



with a more pronounced peak as the limit is ap- 
proached, just as in the double unresolved and 
single soft limits discussed earlier. This demon- 
strates the convergence of the counterterm to the 
matrix element. 

However, the azimuthal correction term 
Q F o(i,j,z,k±) has the unfortunate side effect of 
generating new divergences which are not present 
in the matrix element. 

A second more successful approach is to can- 
cel the angular terms by combining phase space 
points related to each other by a rotation of the 
system of unresolved partons [30113] . When both 
collinear partons are in the final state, this can 
be achieved by considering pairs of phase space 
points which are related by rotating the collinear 
partons by an angle of ir/2 around the resul- 
tant parton direction. Similarly, for the initial- 
final state collinear configurations produced when 
Pi —> +pj for i = 1, 2 and with i\ \j and p 2 = 0, 
the phase space points should again be related by 
rotations of tt/2 about the direction of This 
has the consequence of rotating p^ off the beam 
axis and therefore has to be compensated by a 
Lorentz boost. 

The effect of combining pairs of phase space 
points is shown in Figs. HJc) and EJc). We 
see that the distributions for both final-final and 
initial-final collinear limits have a very sharp 
peak around R = 1. For the final state sin- 
gularity and x = 10~ 12 we obtained an aver- 
age of R = 0.999996 and a standard deviation 
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Figure 1. Distribution of R for 10000 single final 
state collinear phase space points for x = sy/s 
and x — 10~ 8 (red), x = 10 -9 (green) and 
x = 10~ 10 (blue) for (a) the uncorrected matrix 
subtraction term (b) the azimuthally corrected 
subtraction term and (c) pairs of phase space 
points related by a rotation of it/ 2 about the axis 
of the collinear pair. 



Figure 2. Distribution of R for 10000 single 
initial state collinear phase space points for (a) 
the uncorrected matrix subtraction term (b) the 
azimuthally corrected subtraction term and (c) 
pairs of phase space points related by a rotation 
of 7r/2 about the axis of the collinear pair and 
then boosted so the initial direction is preserved. 



of a = 0.00015. Similarly, for the initial state 
collinear singularity and x = — 10~ 12 we found 
R = 0.9999991 and a = 0.00011. There is an 
enormous improvement compared to the raw dis- 



tributions shown in Figs.QJa) and[2](a) and a sig- 
nificant improvement compared to adding an az- 
imuthal correction term shown in Figs, [ljb) and 
[2Jb). The phase space point averaged subtrac- 
tion term clearly converges to the matrix element 
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as we approach the single collinear limit and cor- 
rectly subtracts the azimuthally enhanced terms 
in a point-by-point manner. 

5. Conclusions 

In this contribution, we have discussed the ap- 
plication of the antenna subtraction formalism to 
construct the subtraction term relevant for the 
six-gluon double real radiation contribution to 
di-jet production. The subtraction term is con- 
structed using four-parton and three-parton an- 
tennae. We showed that the subtraction term 
correctly describes the double unresolved limits 
of the gg — > gggg process. In particular, by com- 
bining phase space points, the subtraction term 
avoids the problems associated with angular cor- 
relations produced when a gluon splits into two 
collinear gluons. The final goal is the construction 
of a numerical program to compute NNLO QCD 
estimates of di-jet production in hadron-hadron 
collisions. 
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